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Abstract 

Large-volume, high-temperature Bose-Einstein condensation is illus- 
trated for a relativistic 0(2)-invariant scalar field with fixed charge using 
the canonical ensemble. The standard, grand canonical results are repro- 
duced for the infinite- volume limit. Finite- volume corrections are calcu- 
lated in the canonical ensemble and the results are found to differ from 
the finite-volume grand canonical approximation in a consistent qualita- 
tive way. 

PACS numbers: 03.75Hh, ll.lOWx. 

1 Introduction 

When dealing with conserved charges in a system in equilibrium the usual pro- 
cedure is to employ the grand canonical ensemble 0ElEj- The physical picture 
from which the grand canonical ensemble is derived is that of a system in contact 
with an infinite reservoir of particles with an associated chemical potential for 
flow of particles into and out-of the system. Averaging in the grand canonical 
ensemble is performed over all possible charge states. The average charge is 
then fixed by a definite choice of chemical potential. 

By contrast, in the canonical ensemble, charge cannot move into and out- 
of the system. There is no chemical potential. The charge is fixed and any 
averaging must be performed only over states with this definite fixed charge. 
This provides a more reasonable model for a system which is insulated with 
respect to charge and a better approximation for experiments involving Bose- 
Einstein condensation (BEC) of trapped atoms 0]. 

A clear outline of the path integral formulation of relativistic field theory 
with exactly conserved energy and charge can be found in Refs.0 EJ- The 
central technique is to insert delta-functions into the trace over states in order 
to select only those states with chosen eigenvalues of the Hamiltonian and charge 
operators. This is the technique we shall use to exactly conserve charge in the 
canonical ensemble. An example of the application of this method to BEC in the 
canonical ensemble can be found in Ref. 7 where fluctuations in the ground state 
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occupation are studied in the non-interacting, non-relativistic limit. A further 
study of the comparison between canonical and grand canonical ensembles with 
regard to trapped atoms can be found in Ref. 8 . 

In this article, we outline a fully relativistic calculation of BEC using the 
canonical ensemble. After mathematically formulating the model in the next 
section we go on to consider the high-temperature, large-volume limit in sec- 
tion (c.f. Refs.p |21 El)- We calculate the proportion of the total charge 
that occupies the lowest energy state in section^ and as expected, the canon- 
ical and grand canonical ensembles are found to give the same results in the 
infinite- volume limit. In section we consider finite- volume corrections where 
a difference between the two ensembles is observed. 



2 Exactly conserved charge 

To study a conserved charge Q using the canonical ensemble we must identify 
and only average over those states with this specific charge. As noted above, 
this is taken care of by inserting a delta-function into the trace over states. The 
partition function is then 

Z = Tr K,Q?}- W 

The density operator p in the canonical ensemble gives a probability to each 
state in accordance with the temperature of the system T — 1/(3. It is given by 

p = e-?», (2) 

where H is the Hamiltonian operator. 

The total charge Q has integer value and so the integral form of the delta- 
function is 

(3) 



2vr J_ n 



The partition function now contains a trace over states where the density oper- 
ator is multiplied by a charge operator-dependent phase factor 

Z = J da c lQa Tr je"^ 7 p} . (4) 

This phase factor can be absorbed into the states where its effect is felt through 
a change in the boundary conditions. However, it is simpler to absorb it as a 
shift in the Hamiltonian of the system: 

H->H' = H + i^Q. (5) 

P 

The non-hermiticity of this Hamiltonian would concern us if we were to apply 
it to a dynamical problem. In thermal equilibrium this is nothing more than a 
mathematical trick. 

We shall consider a relativistic (9(2)-invariant scalar theory. This is one 
of the simplest theories which contains a conserved charge. It is useful as a 
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simplified version of the Higgs sector and also for describing atomic gases in its 
non-relativistic limit. The Hamiltonian and charge are given by 



H 



d 3 x- 
v 2 



7T? + ^2 + (V0i) 2 + (V0 2 ) 2 + m 2 ^ + m 2 (j>l 



Q = d 3 X (</> 2 7Tl - 4>l7T 2 ). 



(0) 



The trace in Eq. Q can be expressed as a functional integral in terms of the 
Lagrangian for this system [H] 



(7) 



Tr {e"** p] = Tr {e"^'} oc / [d(/> a ] e^o"* dt L ' (t) 
where the Lagrangian is found to be 

L'(t) = J^x\ 



dt4>i - i-p<t>2 



Qt<i>2 + i-3<t>l 



(V<^) 2 - (V<^ 2 ) 2 - m 2 (fi - m 2 (j>l 



(8) 



To calculate the functional integral we begin by Fourier decomposing the fields 9 . 
Whilst doing this we shall anticipate the existence of a condensate by giving the 
fields a spatially constant component v a to be determined |2J. The constant v a 
acts as a classical replacement for a; o(O). We define r = it where < r < (3. 
The fields are periodic over this imaginary time interval and we have 



Q (x, r ) = „ Q + -L £^ e-» T +lp ' x o; „(p), 
where to n = 2irn//3, and 

d 3 p 

?>)3' 



The action is now given by 

-i/3 



Jo 



dt L'(t) = --I3V m 



P 2 



where, after defining u> 2 = p 2 + m , the inverse propagator is 



- K n (p) 



< + u + ; 

, -2*§^ 2 



2i^ui n 

UJ n + uj 2 + 



(9) 



(10) 



(11) 



(12) 



It can be checked that this action is real. The trace in Eq.Q now becomes 



Tr 



{e-^ p} 



tx e 



tx e 



-i/3y(m 2 +||)|H s 



-i/3v(m 2 + |!)H 2 -ln Vdctif 



32 S l4>a:- n (-p)K ab . n (p)4> b ; n (.p) 



(13) 
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A few lines of algebra determine the integral form for In y/det K. We find 



In VdetK = V [pw + ln(l - e - 0u+ia ) + ln(l - e^" - * 7 )] . (14) 
J P 

The first term on the right-hand side of Eq. (|14Jl is the infinite zero-point con- 
tribution which can be factored out since it doesn't depend on a. The other 
contributions cannot be easily calculated. We are forced to consider their be- 
havior in extreme-temperature limits. Our Ea. (|14f> is essentially the same as 
Eq.(2.57) in Ref.|10| with the replacement a = i/3[i where /i is the chemical 
potential. 



3 High-temperature and large-volume 



Here we shall demonstrate how the contributions to In Vdet K can be calculated 
in the high- temperature and large- volume limit. We first perform a Taylor 
expansion: 



ln(l - e-^ ±,; ") = -V . (15) 

Z — J 77 



The momentum integration can now be taken to act only on the relevant expo- 
nential factor 

/ ln(l - e-^ ±ta ) = -J2 ~ / e "™^- ( 16 ) 

JP n=l U ->V 

To perform this integral in the limit of small [3 we first write out the Mellin- 
Barnes integral form for the exponential 



1 

2ttI 



c+ioo 

dsa- s r(s), (17) 



J ^ = ^ J - ^ {11(3) ' S J 



U> 

P 



The momentum integral is now straightforward 



Ul = 



(4tt)I r(f) 



T(s). (18) 



(19) 



The value of c must be greater than 3 such that the Mellin-Barnes integrand is 
finite. The contour may then be closed at infinity in the negative real half-plane 
where the contribution to the integral is zero. The leading pole comes from the 
factor T (^j^)- Evaluating its residue gives 

/~^~M^) 3 - (20) 

This result could also have been found by taking u) ~ p + • • • . This is a valid 
approximation when [3m <C 1. However, the Mellin transform allows us to 
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consider corrections to this limit in a controlled way (see later). We therefore 
have 

Ml-e-^)~-^£— (21) 

n— 1 



Finally, summing the two logarithmic contributions to In V det K we find |llj 

r ov °° 
V / [ln(l - e -^ +lrT ) + ln(l - e^"-**)] ^— £ ■ 

(22) 



i 4 

r 2„.2 _ I _ 13 ^4 



tt 2 /3 3 V90 12 12 48 

We can now bring everything together in order to write out the canonical par- 
tition function with a conserved charge Q. We may factor out any tr-independent 
terms as an overall normalization TV. We also write in terms of a new variable 
x where 2nx = a. Wc thus have 

Z ~ JV(P) r ^ e -^ y ( m2+ ^) H2 e-^ 2 ^ 2 . (23) 



Within this range of integration, in the limit of large V and small j3 we may 
make the approximation 



_ 21^2 ( | 3,1-1)2 _21^_ x . 

e Vl 1 y — e S S . (24) 

Since this exponential factor suppresses any potential contribution for larger 
values of x we may further extend the integration range to infinity with negligible 
effect giving 

Z ~ Af(p) r dx cos(2^Qx) e -iP v ( m2 +^)^ 2 e -^ x \ (25) 



This is the canonical partition function for an 0(2)-invariant scalar field theory 
with total conserved charge Q. The integral can be performed analytically. In 
the next section we shall use Z to demonstrate Bose-Einstein condensation in 
the canonical ensemble. 



4 Bose-Einstein condensation 

The classical variable v a represents the constant expectation value of the field 
and is to be determined. We first identify \nZ(v a ) with the effective potential 
for the field. The correct solution v a must minimize the effective potential: 

*]nZ=±^=Q. (26) 
dv a Z dv a 

Applying this condition to Z in Ea. l25J) we find 

< , i ;,-d»„e» w ^.-^(^1-- 
1 1 
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This can be evaluated with the following standard integrals 

1 R 



dx cos(ax) e' bx ' = |W~ e~%r, (28) 
d,cosM^e- 2 .i^(l-^)e-*. (29) 
Given these integrals we find 



= v a 



m 2 ; — ] . (30) 

|2 X 



l y2 {w 



This constraint equation looks very much like the grand canonical constraint 
(see Eq.(1.42) in [3]) once we define the chemical potential by 

M = r~ -• (31) 

There are two solutions: 

v a = (32) 

The first solution is trivial. The second solution exists when > ^i^. Wc 
may therefore define a critical temperature to mark the point at which the field 
develops a non-zero expectation value 



T e =±- = J*. (34) 
(3 C V Vm V ' 

To express the charge stored in the ground state we write out the total charge 
in integral form: 



(35) 



giving 

Q = — 5- \Tm+\ v \ — 2V7r 2, i , , , 2 n - (36) 

= Qfluct + Qcond- (37) 

Here we can divide the contributions to the charge into two types. The first term 
in brackets corresponds to the charge contained in the fluctuations or excited 
modes. The second term is the charge which is contained in the ground state, 
the Bosc-Einstcin condensate. We find 

O \v? - f° whe n/ 3 2 <^ 



3/32 when (3 2 > 3Q 
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At T = 1/(3 = all the charge resides in the condensate. As the temperature 
increases charge is excited out. Eventually, when T > T c , the condensate melts 
and all the charge resides in the fluctuation modes. Of course, these equations 
are only valid in the high-temperature limit and we cannot rely on them when 
the condensate becomes big. 

It has been demonstrated how we can understand BEC without using a 
chemical potential. We achieve the same result as for the grand canonical en- 
semble in the high-temperature, large-volume limit. 



5 Finite-volume 



Up until now, the canonical ensemble has proved to be more complicated than 
the grand canonical ensemble without offering any new results. However, now 
we stand to benefit from the fact that the canonical ensemble provides a more 
realistic picture of a finite-volume system, closed to the movement of charge 
into and out-of that volume. To continue from the previous section we work in 
the high-temperature region and the bosons are not subjected to any external 
potential within their volume. 

It is well known that there is no phase transition involved in Bose-Einstein 
condensation for finite volumes. However, we may still characterize the process 
by studying the proportion of the total charge stored in the lowest energy state. 
In the grand canonical ensemble the charge in the lowest energy state is given 



by 



Qo 



2T/x 



(39) 



For T < T c when the charge occupying the lowest energy state compares with 
the charge stored in fluctuations, we may think of (m — fi)^ 1 ~ 0(VmT) \1'2\ . 
The chemical potential is determined by the constraint that the total energy is 
fixed. To leading order in the large-volume expansion we have |13| 



21> 

m 2 — y« 2 



V 



fiT 2 fj,T(m 2 - n 2 )* fi(3m 2 - 2^i 2 ) 



2vr 



12tt 2 



(40) 



For any given T, V and Q, this equation can be solved for /i enabling Qo to be 
calculated. (For a more detailed finite- volume calculation see |14j.) 

Returning to the canonical ensemble and the expression for In \/det K given 
in equation l(Ti|) . for finite volumes, the integral over momentum states should 
be replaced with an infinite sum over those momentum states compatible with 
the boundary conditions. We again expand as follows 



EMi 



n=l 



p 



T10UJ 



1 pc-\-ioo 



r( s ). 



(41) 



(42) 



Choosing a cubic container for our system with side L and Neumann boundary 
conditions we have 
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E 

1 



/7rni\ 2 /7rrii\ 2 /7rni> 2 



\ L J 



L 



dt f2 -i e 



l„-m' 



I 



The sum over p can be performed as a large-volume expansion |14j 



n=0 



(tt) * 



+ 



i 



This gives 



E< 



r 



8(7rf)i Svrf 8(7ri)3 8 

- 3 " sr (^) , T ,|3m^r(^) 



8rrr(§) 



and eventually 



InVdetX V 



87T2T I 



2 cos(ner) , 1 cos(na) 

E— ^ + ^™ 2 ^E 



-V* 



7T 2 /3 3 

3 



n=l 

oo 



2tt 2 /3^ 

n— 1 



Ecos(ner) a 2 o y-^ 
„ 3 



3 

47T 



cos(ner) 



n=l 



-vi—Y 

n—1 



cos(nfj) 



^E 



-n/3m cos ( ncr ) 



(43) 



(44) 



(45) 



(46) 



We have kept all terms in the l/L expansion up to O(L ) and in the high 
temperature (small mass) expansion up to 0(m 2 ). All but one of these infinite 
sums have analytic solutions for the given integration range of a. Referring 
to we have 



cos(ncr) 7r 4 7r 2 er 
n 4 " ~ 90 12~ 



E 



E 

71=1 
OO 

E u- 

n=l 

oo , s 

Ecos(ncr) 
n 

n=l 

e -«/3m cos ( ncr ) 



2 - 2 ttI^I 3 a 4 



12 



48 



cos(ncr) 7r 2 7r|cr| a 1 
" ~6 2~~ + T' 



•-b(2 - 2coscr) 



(47) 
(48) 
(49) 



-- [ln(l - e"' 3m+l<T ) + ln(l - C - f3m - ia )] . (50) 



The remaining sum is convergent and can easily be evaluated numerically. 
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There is no guarantee that the expansion of Eq. l|46|l is a good one. The 
best we can expect is that subsequent terms in the expansion quickly become 
negligibly small in comparison with the leading order terms over the whole 
integration range of a. This is clearly not the case for the fourth term in the 
expansion when close to a = 0. The term diverges leading to a zero in the 
integrand. However, by appropriate choice of parameters, the domination of 
this term is made sufficiently brief that its effect is negligible. 

The O(V ) term of Ea. l|46() may be rewritten as 

\ £ In [(^ + m 2 + T 2 a 2 ) 2 - TW n ] (51) 

n 

from which we may extract the n = contribution and identify this as the 
contribution to In Vdet K of the ground state mode (j> a: Q(0): 

ln(m 2 + TV). (52) 

The expectation of charge in the lowest energy state is given by 

(Qo) = ^f^ clQ ° Tr P Qo} ■ (53) 

This can be derived by differentiating with respect to a only those terms in Z 
where a is coupled to the ground state mode (expression (|52|l ). We find 

where In Vdet K is given above in H46JI • 

Having defined the condensate by the expectation of charge in the lowest 
energy mode, we may directly calculate the condensate fraction by numerical 
integration. The results can be seen for different values of the total charge 
in Fig^ and Fig|2] A comparison is made with the infinite-volume limit and 
with the finite-volume approximation in the grand canonical ensemble given by 
Ea. (|40|l . The parameters must be chosen such that (3 /V 1 ^ 3 , (3m < 1 in order 
to justify the expansion of In \J det K. We chose Q ~ 10 4 , V = 10 9 /T 3 and 
m = 3Q I (VT 2 ) , working in units where T c — 1. 

For T > T c we consistently observe a lower occupation of the ground state 
in the canonical ensemble as compared with the grand canonical ensemble. This 
makes the infinite-volume limit a better approximation for the canonical than 
the grand canonical case. As the temperature approaches the critical tempera- 
ture we see the ground state occupation increase as expected. At low temper- 
atures the expansion of In v det K breaks down and canonical results deviate 
from the other approximations. 

6 Conclusions 

It has been shown how we can understand BEC for a system of fixed charge in 
the canonical ensemble. Though the chemical potential plays a crucial role in 
our understanding of BEC in the grand canonical ensemble, there is no need for 
a chemical potential in the canonical ensemble. 
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Figure 1 : Expectations of charge in the lowest energy state against temperature. The 
solid line represents the infinite-volume approximation in the high-temperature limit; 
the dotted line is a finite- volume approximation using the grand canonical ensemble; 
the dashed line uses the canonical ensemble. Parameters are Q = 20000, V = fO 9 , 
mass is chosen to give T c = 1 (m = 0.00006). 

We have seen that the two ensembles produce identical results in the infinite- 
volume limit. For finite volumes, the canonical method produces sensible results 
which differ from the grand canonical approximation in a consistent qualitative 
way. 

For future work, a low-temperature, non-relativistic version of this calcu- 
lation could feasibly be tested against experimental data from BEC in atomic 
fluids and other theoretical results |S] . 

Acknowledgements 

This work was financially supported by The Royal Commission for the Exhi- 
bition of 1851. I would like to acknowledge Tim Evans and Marko Ivin for 
discussions concerning the canonical ensemble. I would also like to acknowledge 
the hospitality of the particle theory group at the University of Sussex and the 
HEP group at UCL. 

References 

[1] H. E. Haber and H.A. Weldon, Phys. Rev. D25, 502 (1982). 
[2] J. I. Kapusta, Phys. Rev. D24 426 (1981). 

[3] K. M. Benson, J. Bernstein, and S. Dodclson, Phys. Rev. D44 2480 (1991). 



10 




Figure 2: Q = 40000, V = 10 9 , m = 0.00012 in units of T c . Canonical (dashed line), 
grand canonical (dotted line), infinite volume approximation (solid line). 

[4] F. Dalfovo and S. Giorgini, Rev. Mod. Phys. 71 463 (1999). 

[5] M. Chaichian and I. Scnda, Nucl. Phys. B396 737 (1993). 

[6] M. Chaichian, H. Satz, and I. Scnda, Phys. Rev. D49 1566 (1994). 

[7] H. D. Politzcr, Phys. Rev. A54 5048 (1996). 

[8] C. Herzog and M. Olishanii, Phys. Rev. A55 3254 (1997). 

[9] N. P. Landsman and Ch. G. van Weert, Phys. Rep. 145, 141 (1987). 

[10] J. I. Kapusta, Finite- Temperature Field Theory (Cambridge Univ. Press, 
Cambridge 1989). 

[11] I. S. Gradshteyn and I. M. Ryzhik, Tables of integrals, series and products 
(Academic Press 1994). 

[12] J. D. Smith and D. J. Toms, Phys. Rev D53 5771 (1996). 

[13] H. E. Habcr and H. A. Wcldon, Phys. Rev. Lett. 46 1497 (1981). 

[14] K. Shiokawa and B. L. Hu, Phys. Rev. D60 105016 (1999). 



11 



